Innovative direct introduction-ion mobility-mass spectrometry (DI-IM-MS) approach for fast and robust isomer-specific quantification in a complex matrix: Application to 2'-fucosyllactose (2'-FL) in breast milk.

Identification and specific quantification of isomers in a complex biological matrix by mass spectrometry alone is not an easy task due to their identical chemical formula and therefore their same mass-to-charge ratio (m/z). Here, the potential of direct introduction combined with ion mobility-mass spectrometry (DI-IM-MS) for rapid quantification of isomers as human milk oligosaccharides (HMOs) was investigated. Differences in HMO profiles between various analyzed breast milk samples were highlighted using the single ion mobility monitoring (SIM2) acquisition for high ion mobility resolution detection. Furthermore, the Se+ (secretor) or Se- (non-secretor) phenotype could be assigned to breast milk samples studied based on their HMO contents, especially on the response of 2'-fucosyllactose (2'-FL) and lacto-N-fucopentaose I (LNFP I). The possibility of quantifying a specific isomer in breast milk by DI-IM-MS was also investigated. The standard addition method allowed the determination of the 2'-FL despite the presence of other oligosaccharides, including 3-fucosyllactose (3-FL) isomer in breast milk. This proof-of-concept study demonstrated the high potential of such an approach for the rapid and convenient quantification of isomers in complex mixtures.

Assessment of local and systemic signature of eosinophilic esophagitis (EoE) in children through multi-omics approaches.

Background: Eosinophilic oesophagitis (EoE) is a chronic food allergic disorder limited to oesophageal mucosa whose pathogenesis is still only partially understood. Moreover, its diagnosis and follow-up need repeated endoscopies due to absence of non-invasive validated biomarkers. In the present study, we aimed to deeply describe local immunological and molecular components of EoE in well-phenotyped children, and to identify potential circulating EoE-biomarkers. Methods: Blood and oesophageal biopsies were collected simultaneously from French children with EoE (n=17) and from control subjects (n=15). Untargeted transcriptomics analysis was performed on mRNA extracted from biopsies using microarrays. In parallel, we performed a comprehensive analysis of immune components on both cellular and soluble extracts obtained from both biopsies and blood, using flow cytometry. Finally, we performed non-targeted plasma metabolomics using liquid chromatography coupled to high-resolution mass spectrometry (LC-HRMS). Uni/multivariate supervised and non-supervised statistical analyses were then conducted to identify significant and discriminant components associated with EoE within local and/or systemic transcriptomics, immunologic and metabolomics datasets. As a proof of concept, we conducted multi-omics data integration to identify a plasmatic signature of EoE. Results: French children with EoE shared the same transcriptomic signature as US patients. Network visualization of differentially expressed (DE) genes highlighted the major dysregulation of innate and adaptive immune processes, but also of pathways involved in epithelial cells and barrier functions, and in perception of chemical stimuli. Immune analysis of biopsies highlighted EoE is associated with dysregulation of both type (T) 1, T2 and T3 innate and adaptive immunity, in a highly inflammatory milieu. Although an immune signature of EoE was found in blood, untargeted metabolomics more efficiently discriminated children with EoE from control subjects, with dysregulation of vitamin B6 and various amino acids metabolisms. Multi-blocks integration suggested that an EoE plasma signature may be identified by combining metabolomics and cytokines datasets. Conclusions: Our study strengthens the evidence that EoE results from alterations of the oesophageal epithelium associated with altered immune responses far beyond a simplistic T2 dysregulation. As a proof of concept, combining metabolomics and cytokines data may provide a set of potential plasma biomarkers for EoE diagnosis, which needs to be confirmed on a larger and independent cohort.

Rapid structural characterization of human milk oligosaccharides and distinction of their isomers using trapped ion mobility spectrometry time-of-flight mass spectrometry.

Oligosaccharides have multiple functions essential for health. Derived from the condensation of two to several monosaccharides, they are structurally diverse with many co-occurring structural isomer families, which make their characterization difficult. Thanks to its ability to separate small molecules based on their mass, size, shape, and charge, ion mobility-mass spectrometry (IM-MS) has emerged as a powerful tool for separating glycan isomers. Here, the potential of such a technique for the rapid characterization of main human milk oligosaccharides (HMOs) was investigated. Our study focused on 18 HMO standards. The IM-MS analysis enabled to distinguish almost all the HMOs studied, in particular thanks to the single ion mobility monitoring acquisition using the trapped ion mobility spectrometry device, providing high ion mobility resolution and enhanced ion mobility separation. Alternatively, the combination of IM-MS separation with MS/MS experiments has proven to increase performance in identifying HMOs and especially isomers poorly separated by ion mobility alone. Finally, collision cross-section values are provided for each species generated from the 18 HMOs standards, which can serve as an additional identifier to characterize HMOs.

Long-range linkage disequilibrium in French beef cattle breeds.

BACKGROUND: Linkage disequilibrium (LD) is a key parameter to study the history of populations and to identify and fine map quantitative trait loci (QTL) and it has been studied for many years in animal populations. The advent of new genotyping technologies has allowed whole-genome LD studies in most cattle populations. However, to date, long-range LD (LRLD) between distant variants on the genome has not been investigated in detail in cattle. Here, we present the first comprehensive study of LRLD in French beef cattle by analysing data on 672 Charolais (CHA), 462 Limousine (LIM) and 326 Blonde d'Aquitaine (BLA) individuals that were genotyped on the Illumina BovineHD Beadchip. Furthermore, whole-genome LD and haplotype block structure were analysed in these three breeds. RESULTS: We computed linkage disequilibrium (r2) values for 5.9, 5.6 and 6.0 billion pairs of SNPs on the 29 autosomes of CHA, LIM and BLA, respectively. Mean r2 values drop to less than 0.1 for distances between SNPs greater than 120 kb. However, for the first time, we detected the existence of LRLD in the three main French beef breeds. In total, 598, 266, and 795 LRLD events (r2 ≥ 0.6) were detected in CHA, LIM and BLA, respectively. Each breed had predominantly population-specific LRLD interactions, although shared LRLD events occurred in a number of regions (55 LRLD events were shared between two breeds and nine between the three breeds). Examples of possible functional gene interactions and QTL co-location were observed with some of these LRLD events, which suggests epistatic selection. CONCLUSIONS: We identified long-range linkage disequilibrium for the first time in French beef cattle populations. Epistatic selection may be the main source of the observed LRLD events, but other forces may also be involved. LRLD information should be accounted for in genome-wide association studies.

Gene networks for three feed efficiency criteria reveal shared and specific biological processes.

BACKGROUND: French beef producers suffer from the decrease in profitability of their farms mainly because of the continuous increase in feed costs. Selection for feed efficiency in beef cattle represents a relevant solution to face this problem. However, feed efficiency is a complex trait that can be assessed by three major criteria: residual feed intake (RFI), residual gain (RG) and feed efficiency ratio (FE), which involve different genetic determinisms. An analysis that combines phenotype and whole-genome sequence data provides a unique framework for genomic studies. The aim of our study was to identify the gene networks and the biological processes that are responsible for the genetic determinism that is shared between these three feed efficiency criteria. RESULTS: A population of 1477 French Charolais young bulls was phenotyped for feed intake (FI), average daily gain (ADG) and final weight (FW) to estimate RFI, RG and FE. A subset of 789 young bulls was genotyped on the BovineSNP50 single nucleotide polymorphism (SNP) array and imputed at the sequence level using RUN6 of the 1000 Bull Genomes Project. We conducted a genome-wide association study (GWAS) to estimate the individual effect of 8.5 million SNPs and applied an association weight matrix (AWM) approach to analyse the results, one for each feed efficiency criterion. The results highlighted co-association networks including 626 genes for RFI, 426 for RG and 564 for FE. Enrichment assessment revealed the biological processes that show the strongest association with RFI, RG and FE, i.e. digestive tract (salivary, gastric and mucin secretion) and metabolic processes (cellular and cardiovascular). Energetic functions were more associated with RFI and FE and cardio-vascular and cellular processes with RG. Several hormones such as apelin, glucagon, insulin, aldosterone, the gonadotrophin releasing hormone and the thyroid hormone were also identified, and these should be tested in future studies as candidate biomarkers for feed efficiency. CONCLUSIONS: The combination of network and pathway analyses at the sequence level led to the identification of both common and specific mechanisms that are involved in RFI, RG and FE, and to a better understanding of the genetic determinism underlying these three criteria. The effects of the genes involved in each of the identified processes need to be tested in genomic evaluations to confirm the potential gain in reliability of using functional variants to select animals for feed efficiency.

Progressive ataxia of Charolais cattle highlights a role of KIF1C in sustainable myelination.

Hereditary spastic paraplegias (HSPs) are clinically and genetically heterogeneous human neurodegenerative diseases. Amongst the identified genetic causes, mutations in genes encoding motor proteins such as kinesins have been involved in various HSP clinical isoforms. Mutations in KIF1C are responsible for autosomal recessive spastic paraplegia type 58 (SPG58) and spastic ataxia 2 (SPAX2). Bovines also develop neurodegenerative diseases, some of them having a genetic aetiology. Bovine progressive ataxia was first described in the Charolais breed in the early 1970s in England and further cases in this breed were subsequently reported worldwide. We can now report that progressive ataxia of Charolais cattle results from a homozygous single nucleotide polymorphism in the coding region of the KIF1C gene. In this study, we show that the mutation at the heterozygous state is associated with a better score for muscular development, explaining its balancing selection for several decades, and the resulting high frequency (13%) of the allele in the French Charolais breed. We demonstrate that the KIF1C bovine mutation leads to a functional knock-out, therefore mimicking mutations in humans affected by SPG58/SPAX2. The functional consequences of KIF1C loss of function in cattle were also histologically reevaluated. We showed by an immunochemistry approach that demyelinating plaques were due to altered oligodendrocyte membrane protrusion, and we highlight an abnormal accumulation of actin in the core of demyelinating plaques, which is normally concentrated at the leading edge of oligodendrocytes during axon wrapping. We also observed that the lesions were associated with abnormal extension of paranodal sections. Moreover, this model highlights the role of KIF1C protein in preserving the structural integrity and function of myelin, since the clinical signs and lesions arise in young-adult Charolais cattle. Finally, this model provides useful information for SPG58/SPAX2 disease and other demyelinating lesions.

Insights into the genetic variation of maternal behavior and suckling performance of continental beef cows.

BACKGROUND: In beef cattle, maternal care is critical for calf survival and growth. Our objective was to evaluate the major sources of additive genetic variation in maternal behavior and suckling performance in two genetically close beef breeds. METHODS: Maternal performance was assessed based on maternal behavior (MB), milk yield (MY) and udder swelling score (US) of 1236 Blonde d'Aquitaine cows and 1048 Limousin cows. MB was scored just after calving to describe the intensity of the dam's protective behavior towards her calf. Most of the cows were genotyped using the low-density chip EuroG10K BeadChip, and imputed to the high-density 770K panel within breed. Genetic parameters for each trait were estimated for each breed under a multi-trait best linear unbiased prediction animal model. Genomic analysis was performed for each breed using the high-density genotypes and a Bayesian variable selection method. RESULTS: Heritabilities were low for MB (0.11-0.13), intermediate for MY (0.33-0.45) and high for US (0.47-0.64). Genetic correlations between the traits ranged from 0.31 to 0.58 and 0.72 to 0.99 for the Blonde d'Aquitaine and Limousin breeds, respectively. Two quantitative trait loci (QTL) were detected for MB in Blonde d'Aquitaine with NPY1R and ADRA2A as candidate causative genes. Thirty to 56 QTL were detected for MY and US in both breeds and 12 candidate genes were identified as having a role in the genetic variation of suckling performance. Since very few pleiotropic QTL were detected, there was little biological explanation for the moderate (0.57) to very high (0.99) genetic correlations estimated between MY and US in the Blonde d'Aquitaine and Limousin cows, respectively. In Blonde d'Aquitaine, the correlation was largely due to the pleiotropic QTL detected in the region upstream of the CG gene, while in Limousin, this region was only identified for US, thus attesting the difference in genetic architecture between the breeds. CONCLUSIONS: Our findings question the assumption that two populations that have close genetic links share many QTL. Nevertheless, we identified four candidate genes that may explain a substantial amount of the genetic variation in suckling performance of these two breeds.

Detection of quantitative trait loci for maternal traits using high-density genotypes of Blonde d'Aquitaine beef cattle.

BACKGROUND: The genetic determinism of the calving and suckling performance of beef cows is little known whereas these maternal traits are of major economic importance in beef cattle production systems. This paper aims to identify QTL regions and candidate genes that affect maternal performance traits in the Blonde d'Aquitaine breed. Three calving performance traits were studied: the maternal effect on calving score from field data, the calving score and pelvic opening recorded in station for primiparous cows. Three other traits related to suckling performance were also analysed: the maternal effect on weaning weight from field data, milk yield and the udder swelling score recorded in station for primiparous cows. A total of 2,505 animals were genotyped from various chip densities and imputed in high density chips for 706,791 SNP. The number of genotyped animals with phenotypes ranged from 1,151 to 2,284, depending on the trait considered. RESULTS: QTL detections were performed using a Bayes C approach. Evidence for a QTL was based on Bayes Factor values. Putative candidate genes were proposed for the QTL with major evidence for one of the six traits and for the QTL shared by at least two of the three traits underlying either calving or suckling performance. Nine candidate genes were proposed for calving performance among the nine highlighted QTL regions. The neuroregulin gene on chromosome 27 was notably identified as a very likely candidate gene for maternal calving performance. As for suckling abilities, seven candidate genes were identified among the 15 highlighted QTL. In particular, the Group-Specific Component gene on chromosome 6, which encodes vitamin D binding protein, is likely to have a major effect on maternal weaning weight in the Blonde d'Aquitaine breed. This gene had already been linked to milk production and clinical mastitis in dairy cattle. CONCLUSION: In the near future, these QTL findings and the preliminary proposals of candidate genes which act on the maternal performance of beef cows should help to identify putative causal mutations based on sequence data from different cattle breeds.

High-density marker imputation accuracy in sixteen French cattle breeds.

BACKGROUND: Genotyping with the medium-density Bovine SNP50 BeadChip® (50K) is now standard in cattle. The high-density BovineHD BeadChip®, which contains 777,609 single nucleotide polymorphisms (SNPs), was developed in 2010. Increasing marker density increases the level of linkage disequilibrium between quantitative trait loci (QTL) and SNPs and the accuracy of QTL localization and genomic selection. However, re-genotyping all animals with the high-density chip is not economically feasible. An alternative strategy is to genotype part of the animals with the high-density chip and to impute high-density genotypes for animals already genotyped with the 50K chip. Thus, it is necessary to investigate the error rate when imputing from the 50K to the high-density chip. METHODS: Five thousand one hundred and fifty three animals from 16 breeds (89 to 788 per breed) were genotyped with the high-density chip. Imputation error rates from the 50K to the high-density chip were computed for each breed with a validation set that included the 20% youngest animals. Marker genotypes were masked for animals in the validation population in order to mimic 50K genotypes. Imputation was carried out using the Beagle 3.3.0 software. RESULTS: Mean allele imputation error rates ranged from 0.31% to 2.41% depending on the breed. In total, 1980 SNPs had high imputation error rates in several breeds, which is probably due to genome assembly errors, and we recommend to discard these in future studies. Differences in imputation accuracy between breeds were related to the high-density-genotyped sample size and to the genetic relationship between reference and validation populations, whereas differences in effective population size and level of linkage disequilibrium showed limited effects. Accordingly, imputation accuracy was higher in breeds with large populations and in dairy breeds than in beef breeds. More than 99% of the alleles were correctly imputed if more than 300 animals were genotyped at high-density. No improvement was observed when multi-breed imputation was performed. CONCLUSION: In all breeds, imputation accuracy was higher than 97%, which indicates that imputation to the high-density chip was accurate. Imputation accuracy depends mainly on the size of the reference population and the relationship between reference and target populations.

Novel insights into the bovine polled phenotype and horn ontogenesis in Bovidae.

Despite massive research efforts, the molecular etiology of bovine polledness and the developmental pathways involved in horn ontogenesis are still poorly understood. In a recent article, we provided evidence for the existence of at least two different alleles at the Polled locus and identified candidate mutations for each of them. None of these mutations was located in known coding or regulatory regions, thus adding to the complexity of understanding the molecular basis of polledness. We confirm previous results here and exhaustively identify the causative mutation for the Celtic allele (PC) and four candidate mutations for the Friesian allele (PF). We describe a previously unreported eyelash-and-eyelid phenotype associated with regular polledness, and present unique histological and gene expression data on bovine horn bud differentiation in fetuses affected by three different horn defect syndromes, as well as in wild-type controls. We propose the ectopic expression of a lincRNA in PC/p horn buds as a probable cause of horn bud agenesis. In addition, we provide evidence for an involvement of OLIG2, FOXL2 and RXFP2 in horn bud differentiation, and draw a first link between bovine, ovine and caprine Polled loci. Our results represent a first and important step in understanding the genetic pathways and key process involved in horn bud differentiation in Bovidae.

A 3.7 Mb deletion encompassing ZEB2 causes a novel polled and multisystemic syndrome in the progeny of a somatic mosaic bull.

Polled and Multisystemic Syndrome (PMS) is a novel developmental disorder occurring in the progeny of a single bull. Its clinical spectrum includes polledness (complete agenesis of horns), facial dysmorphism, growth delay, chronic diarrhea, premature ovarian failure, and variable neurological and cardiac anomalies. PMS is also characterized by a deviation of the sex-ratio, suggesting male lethality during pregnancy. Using Mendelian error mapping and whole-genome sequencing, we identified a 3.7 Mb deletion on the paternal bovine chromosome 2 encompassing ARHGAP15, GTDC1 and ZEB2 genes. We then produced control and affected 90-day old fetuses to characterize this syndrome by histological and expression analyses. Compared to wild type individuals, affected animals showed a decreased expression of the three deleted genes. Based on a comparison with human Mowat-Wilson syndrome, we suggest that deletion of ZEB2, is responsible for most of the effects of the mutation. Finally sperm-FISH, embryo genotyping and analysis of reproduction records confirmed somatic mosaicism in the founder bull and male-specific lethality during the first third of gestation. In conclusion, we identified a novel locus involved in bovid horn ontogenesis and suggest that epithelial-to-mesenchymal transition plays a critical role in horn bud differentiation. We also provide new insights into the pathogenicity of ZEB2 loss of heterozygosity in bovine and humans and describe the first case of male-specific lethality associated with an autosomal locus in a non-murine mammalian species. This result sets PMS as a unique model to study sex-specific gene expression/regulation.